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Abstract 

A swarm  intelligence  based  method  is  presented  in  this  paper  to  solve  the  problem  of  parameterization  of  given 
set  of  planar  and  spatial  data  points,  which  may  be  uniformly  or  non-uniformly  spaced,  for  least  squares  fitting 
with  Bezier  curves.  The  basic  framework  of  the  proposed  method  relies  on  fireworks  algorithm.  Performance  of 
the  proposed  method  is  tested  through  several  example  curves  of  varying  complexities. 
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I.  Introduction 

Fitting  of  parametric  curves  and  surfaces  to  measured  or  generated  point  clouds  finds  applications  in 
many  industrial  and  technological  domains,  such  as  computer-aided  design  and  manufacturing  (CAD/ 
CAM),  medical  imaging,  reverse  engineering,  virtual  reality,  etc.  While  the  point  clouds  are  in  three- 
dimensional  Cartesian  space , the  curves  and  surfaces  are  created  in  parametric  space.  The  process  of 
assigning  unique  parameter  values  for  each  point  in  the  point  cloud  is  known  as  parameterization. 
The  accuracy  of  fitted  curves  and  surfaces  depend  on  proper  allocation  of  parameter  values.  Among 
the  available  parameterization  methods,  uniform,  chord  length  and  centripetal  methods  are  commonly 
used.  Fitting  of  parametric  Bezier  curves  for  the  planar  and  spatial  points  is  considered  in  the  present 
work. 

Bezier  curve  fitting  is  generally  based  on  least  squares  method  (LSM).  Curve  fitting  starts  with  the 
estimation  of  parameters  and  knot  vectors  followed  by  the  determination  of  its  control  points  by 
minimizing  the  sum  of  squares  of  distances  between  data  points  and  the  fitted  curve.  An  attempt  to  fit 
B -spline  curves  for  randomly  measured  points  using  the  least  squares  method  has  been  reported  [1].  It 
may  be  noted  that  the  number  of  unknown  variables  to  be  solved  increases  with  large  numbers  of  data 
points.  In  such  cases,  use  of  conventional  optimization  methods  for  solving  the  system  of  equations 
may  not  be  feasible  or  may  result  in  inferior  curves.  This  has  led  researchers  to  explore  the  possibility 
of  solving  curve  fitting  problems  using  neural  networks  [2],  evolutionary  algorithms , etc.  Some  of  the 
research  efforts  related  to  evolutionary  algorithms  are  presented  below. 

Genetic  algorithm  (GA)  is  a commonly  used  evolutionary  algorithm  for  solving  complex  optimization 
problems.  Kumar,  et  al.  proposed  a parameterization  method  for  fitting  non-uniform  B-spline  curves 
using  genetic  algorithms  and  radial  basis  function  neural  networks  [3].  Galvez,  et  al.  developed  two 
artificial  intelligence  techniques  for  the  least-squares  approximation  of  Bezier  curves  and  surfaces, 
viz.  a genetic  algorithm  based  approach  for  parameterization  and  a functional  networks  scheme  for 
handling  additional  functional  constraints  used  [4].  Sun,  et  al.  proposed  an  adaptive  genetic  algorithm 
to  optimize  the  parameters  of  B-spline  curves,  as  the  randomly  initialized  regular  genetic  algorithm 
typically  lead  to  local  minima  and  require  more  control  points  to  assure  higher  fitting  accuracy  [5]. 
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Zhao,  et  al.  applied  genetic  algorithms  to  optimize  the  parameters  of  Bernstein  basis  functions  for 
least  squares  fitting  of  Bezier  curves  to  the  planar  data  points  [6].  Hasegawa,  et  al.  proposed  a multi- 
objective evolutionary  algorithm,  known  as  generalized  differential  evolution  algorithm  3 (GDE3), 
for  parameter  optimization  in  fitting  of  Bezier  curves.  The  authors  have  used  two  objective  functions, 
viz.  error  of  fit  and  length  of  the  curve,  to  obtain  aesthetically  pleasing  curves  while  achieving  a good 
trade-off  between  the  error  and  length  [7].  Galvez  and  Iglesias  presented  an  iterative  mutually  coupled 
hybrid  evolutionary  approach  for  fitting  of  B-spline  curves.  Their  approach  used  genetic  algorithm 
for  data  parameterization  and  the  particle  swarm  optimization  (PSO)  for  knot  placement  [8]. 

Loucera,  et  al.  attempted  to  solve  the  problem  of  approximating  given  set  of  noisy  and  noise-free  data 
points  with  Bezier  curves  through  a hybrid  strategy  that  combines  linear  least-squares  minimization , 
simulated  annealing  and  information  science  metrics , namely  data-fidelity  and  model-complexity  [9]. 
Galvez,  et  al.  used  the  clonal  selection  theory , an  artificial  immune  system  (AIS)  algorithm , to  solve 
the  problem  of  obtaining  the  rational  Bezier  curve  that  fits  a given  set  of  data  points  better  in  the  least- 
squares  sense  [10].  Iglesias,  et  al.  presented  a hybrid  mesh  adaptive  search  algorithm  that  uses  clonal 
selection  theory  as  first  step,  followed  by  an  efficient  local  search  procedure  [11].  This  approach  has 
been  shown  to  yield  better  results  than  just  using  the  clonal  selection  theory. 

Recently,  the  swarm  intelligence  algorithms  are  being  used  to  obtain  suitable  parameters  for  accurate 
fitting  of  curves  and  surfaces.  Swarm  intelligence  is  the  study  of  computational  systems  inspired  by 
collective  intelligence  [12].  Galvez  and  Iglesias  applied  a powerful  metaheuristic  algorithm , called  the 
firefly  algorithm  inspired  from  the  social  flashing  behaviour  of  fireflies  in  nature,  to  obtain  an  optimal 
approximating  Bezier  curve  to  given  set  of  data  points  with  proper  selection  of  control  parameters 
[13].  The  authors  also  used  the  firefly  algorithm  to  approximate  a given  set  of  noisy  data  points  using 
explicit  B-spline  curves  [14].  Iglesias  and  Galvez  presented  an  approach  to  fit  smooth  rational  Bezier 
curves  for  noisy  data  points  using  memetic  firefly  algorithms  by  formulating  the  fitting  problem  to  be 
a continuous  multivariate  nonlinear  optimization  problem  [15].  Iglesias,  et  al.  presented  an  approach 
based  on  the  bat  algorithm  for  optimal  parameterization  for  least-squares  fitting  of  polynomial  Bezier 
curves  [16].  The  application  of  bat  algorithm  has  also  been  extended  to  rational  Bezier  curves  [17]. 
Iglesias  and  Galvez  used  the  memetic  electromagnetism  algorithm  for  approximating  the  data  points 
with  rational  Bernstein  polynomial  curves.  This  algorithm  combines  a global  optimization  algorithm, 
viz.  electromagnetism  algorithm,  with  a local  search  method  [18]. 

Every  evolutionary  algorithm-based  approach  involves  the  tuning  of  several  parameters  that  form  the 
basic  framework  in  order  to  achieve  good  performance  of  the  particular  algorithm.  The  values  of  these 
parameters  are  usually  set  by  user-experience  and  for  a completely  new  problem  they  are  set  by  large 
collection  of  empirical  results.  Hence,  an  algorithm  which  involves  parameters  that  are  easy  to  set  is 
needed.  One  such  algorithm  is  the  fireworks  algorithm  [19]  that  involves  only  four  parameters.  These 
parameters  are  almost  same  for  every  parameterization  problem.  Some  of  the  essential  details  of  this 
algorithm  are  given  in  Section  III.  In  the  present  work,  a fireworks  algorithm  based  parameterization 
approach  is  proposed  for  fitting  of  polynomial  Bezier  curves  and  its  performance  has  been  tested  with 
several  example  cases. 

Remainder  of  the  paper  is  organized  as  follows.  Section  II  describes  about  the  least  squares  fitting  of 
Bezier  curves  and  Section  III  describes  about  the  proposed  method.  Section  IV  presents  and  discusses 
the  results  obtained  using  the  proposed  method  and  finally  the  conclusions  are  presented  in  Section  V. 

II.  Least  Squares  Fitting  of  Bezier  Curves 

The  Bezier  curve  can  be  mathematically  defined  as  [20]: 

p(u)  = 'fplBlJu)\  0<m<1  (1) 

/=0 

where,  the  vectors  pt  represent  the  (n  + 1)  vertices  of  a characteristic  or  control  polygon  and  Bi,n(u)  is 
the  fh  Bernstein  basis  function  of  degree  n.  The  vertices  are  known  as  control  points  as  they  control 
the  shape  of  the  Bezier  curve.  Eq.  (1)  can  be  used  to  compute  the  points  on  a Bezier  curve.  The  values 
of  Bernstein  basis  function  can  be  computed  using  Eq.  (2). 
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f (1  -uT1 
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Let  us  consider  the  reverse  engineering  scenario,  where  the  data  points  are  known  and  the  control 
points  are  to  be  estimated.  Assume  that  there  are  N data  points  ( Dj ) in  a ^-dimensional  space  (£  = 2 for 
planar  curves  and  3 for  space  curves).  The  control  points  of  the  Bezier  curve  can  be  estimated  using 
the  least  squares  ( LS ) fitting  method,  i.e.  by  minimizing  the  sum  of  squared  errors  ( E ) over  the  set  of 
data  points.  The  LS  fitting  can  be  represented  as  follows: 

E = YJ[DJ-p(u )]2  (3) 

7=1 

The  system  of  equations  to  be  solved  by  LS  fitting  can  be  represented  using  Eq.  (4). 

D.i  = YjQ,B1Ju%  0<m<1;  j = 0,...lN-l)  (4) 

1=0 

The  above  equation  can  be  written  in  matrix  form  as  shown  below. 

[D]  = [B][Q]  (5) 

In  Eq.  (5),  vector  [D]  contains  the  data  points,  vector  [ Q ] contains  the  control  points  and  matrix  [B] 
contains  the  Bernstein  basis  function  matrix.  The  control  points  [Q\  may  be  obtained  using  Eq.  (6). 

[Q\  = [BTl[D]  (6) 


III.  Proposed  Method 

In  solving  the  system  of  equation  in  Eq.  (4)  and  Eq.  (10),  suitable  values  of  u (for  curves)  and  u and  v 
(for  surfaces)  are  to  be  determined  first.  Given  the  nonlinear  nature  of  the  problem,  the  parameters  are 
to  be  determined  through  optimization.  The  fireworks  algorithm  is  used  for  solving  this  system  of 
equations.  The  flowchart  of  the  proposed  method  is  shown  in  Figure  1.  Some  of  the  preparatory  steps 
needed  for  the  fireworks  algorithm  are  as  follows: 

1)  The  parameter  values  are  considered  as  fireworks  and  are  encoded  as  real  coded  vectors  of  length 

equal  to  number  of  given  data  points  (N),  i.e.  u=  [ui,  U2,  m\. 

2)  The  objective  function  is  the  function  E given  in  Eq.  (3). 

3)  The  degree  of  the  underlying  curve  and  number  of  control  points  are  not  known  apriori  and  are 
dependent  on  the  complexity  of  the  curve.  The  method  therefore  starts  with  less  number  of  control 
points,  and  increases  it  until  the  error  reduces  below  a certain  threshold  value. 

4)  Setup  the  control  parameters,  namely  the  population  size  ( npop ),  the  number  of  iterations  (n^),  total 
number  of  sparks  (nspr)  and  the  maximum  explosion  amplitude  (A). 


3.1  Number  of  sparks 

The  number  of  sparks  to  be  generated  for  each  member  in  the  population  can  be  found  using  Eq.  (7). 


(7) 


where,  £ is  a small  constant  introduced  to  avoid  the  zero  division  error. 


3.2  Amplitude  of  explosion 

The  amplitude  of  explosion  of  each  member  in  the  population  can  be  estimated  using  Eq.  (8). 


A =A 


^-Emi  n+^r 


(8) 


where,  £ is  a small  constant  introduced  to  avoid  the  zero  division  error. 
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Figure  1.  Flowchart  of  the  proposed  method 


3.3  Generation  of  sparks 

Sparks  are  generated  by  mimicking  the  explosion  process.  Two  kinds  of  spark  generation  methods  are 
generally  employed  considering  the  convergence  and  diversity  of  explosion.  One  method  generates 
the  sparks  towards  the  best  among  the  population,  while  the  other  generates  sparks  in  a random  way. 
Both  ways  of  generating  sparks  are  given  below  in  Eq.  (9). 

u-  = u{  +^.  rand  (Of)  sign(u{best-u{) 

u?  = u{  + (-1)  Wi'(2’U)  ^ rand  (Of)  (9) 

uiGu;  i = 1,2,...  JV; 

where,  the  superscripts  a and /stand  for  sparks  and  fireworks  respectively. 

3.4  Selection  of  sparks  or  locations 

At  the  end  of  every  iteration,  the  current  best  member,  for  which  the  objective  function  is  optimal 
among  population  is  kept  for  the  next  explosion  generation.  After  that,  (npop  - 1)  members  are  selected 
based  on  their  distance  to  other  members  so  as  to  keep  diversity. 

IV.  Results  and  Discussion 

The  chord  length  parameterization  and  proposed  fireworks  algorithm  based  parameterization  methods 
have  been  implemented  using  MATLAB.  The  former  has  been  implemented  for  comparison  purposes. 
Several  example  Bezier  curves  of  varying  complexities  have  been  used  to  compare  these  methods. 
Four  examples  of  planar  curves  and  two  examples  of  space  curves  approximated  by  both  the  methods 
are  shown  respectively  in  Tables  1 and  2.  Axes  of  the  plots  are  marked  in  normalized  coordinates.  A 
visual  examination  of  the  fitted  curves  clearly  reveal  that  the  proposed  parameterization  method  yields 
better  looking  curves  than  those  obtained  using  chord  length  method. 
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Table  1.  Planar  curves  approximated  by  chord  length  and  proposed  parameterization  methods 

Example  Approximation  Method 


No. Chord  Length  Parameterization Proposed  Parameterization 


In  tables  1 and  2,  the  points  marked  with  circle  are  the  data  points  and  those  marked  with  plus  sign  are 
points  on  the  approximated  Bezier  curves.  Table  3 shows  the  errors  associated  with  the  fitted  curves 
using  both  the  methods.  The  dashes  in  Ez  column  for  the  first  four  examples  mean  that  they  are  planar 
curves.  Table  3 clearly  reveals  that  the  proposed  parameterization  method  yields  far  smaller  error-of- 
fit  values  than  those  obtained  using  the  chord  length  parameterization  method.  Thus,  we  can  conclude 
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that  the  proposed  parameterization  method  performs  well  in  terms  of  both  fitting  of  visually  appealing 
Bezier  curves  and  yielding  much  smaller  error-of-fit  values. 

Table  2.  Space  curves  approximated  by  chord  length  and  proposed  parameterization  methods 

Example  Approximation  Method 

No.  Chord  Length  Parameterization  Proposed  Parameterization 


Table  3.  Error-of-fit  values  yielded  by  both  parameterization  methods  for  the  sample  curves  considered 


Example 

Chord  length  parameterization  method 

Proposed  parameterization  method 

No. 

Ex 

Ey 

Ez 

Ex 

Ey 

Ez 

PI 

0.0200 

0.1089 

- 

7.12E-04 

2.36E-04 

- 

P2 

0.0050 

0.0648 

- 

6.93E-05 

1.54E-04 

- 

P3 

0.0010 

0.0400 

- 

8.00E-04 

5.00E-05 

- 

P4 

0.0020 

0.0100 

- 

2.80E-04 

4.90E-04 

- 

SI 

0.0012 

0.0033 

8.84E-04 

2.77E-04 

3.64E-04 

4.78E-04 

S2 

0.0037 

0.0040 

5.40E-06 

3.35E-04 

8.40E-04 

8.91E-06 

V.  Conclusions  and  Future  Scope 

Parameterization  is  a crucial  step  in  fitting  of  parametric  curves  and  surfaces.  This  paper  introduces  a 
new  method  to  address  the  parameterization  problem  in  fitting  of  Bezier  curves  by  approximating  the 
given  point  cloud  data  in  the  least  squares  sense.  Based  on  the  experimentation  using  the  point  clouds 
of  curves  of  varying  complexities,  it  has  been  found  that  the  proposed  parameterization  method  yields 
visually  appealing  curves  along  with  very  small  fitting  errors. 

An  extension  of  the  proposed  method  to  fitting  of  the  Bezier  surfaces  is  being  reported  separately. 
Efforts  are  underway  to  extend  the  proposed  parameterization  method  for  fitting  of  B -Spline  curves 
and  surfaces  also. 
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